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Abstract. The evolution speed in projective Hilbert space is considered for 
Hermitian Hamiltonians and for non-Hermitian (NH) ones. Based on the Hilbert- 
Schmidt norm and the spectral norm of a Hamiltonian, resource-related upper 
bounds on the evolution speed are constructed. These bounds are valid also for 
NH Hamiltonians and they are illustrated for an optical NH Hamiltonian and for 
a non-Hermitian VT— symmetric matrix Hamiltonian. Furthermore, the concept 
of quantum speed efficiency is introduced as measure of the system resources 
directly spent on the motion in the projective Hilbert space. A recipe for the 
construction of time-dependent Hamiltonians which ensure 100% speed efficiency 
is given. Generally these efficient Hamiltonians are NH but there is a Hermitian 
efficient Hamiltonian as well. Finally, the extremal case of a non-Hermitian non- 
diagonalizable Hamiltonian with vanishing energy difference is shown to produce 
a 100% efficient evolution with minimal resources consumption. 



1. Introduction 

Non-Hermitian models naturally emerge in many fields of physics as efficient tools for 
the description of complicated large systems in terms of smaller effective subsystems 
[TJ [5]. Examples range from atomic/molecular physics [31 0j, light propagation in 
optically active crystals [5] and media with anisotropic pumping and absorption [HI E1IH1 
[3J [TUl El HH EH HI EE US] over microwave cavities [T71 [T5] , coupled electronic circuits 
|19[ [20] up to mechanics [211 122) . hydrodynamics [23, 24J and magnetohydrodynamics 
[25, 26, 27, 28J. Apart from the spectral properties of the Hamiltonians, the evolution 
processes generated by non-Hermitian (NH) Hamiltonians can significantly differ from 
those generated by Hermitian Hamiltonians [5, 29, 30, 31, 32] . In this regard it appears 
natural to ask what new possibilities NH evolution entails and what bounds can be 
broken when a Hamiltonian is no longer Hermitian. 

For example, it was shown in |29j that a NH 2x2 matrix Hamiltonian with some 
predefined energy difference can generate a much faster evolution than a Hermitian 
Hamiltonian with the same energy difference. The evolution speed is the rate in 
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which a state changes into other states (e.g. the angular speed of the state vector on 
a Bloch sphere). In the Hermitian case, the energy difference sets an upper bound 
on the evolution speed (the Fleming bound [33 ). The implicitly demonstrated clear 
violation of this bound by non-unitary evolution processes of NH systems leads to 
the conclusion that these energy-difference based bounds should be replaced by some 
more adequate bounds for NH systems. 

The first goal of this article is to derive an upper bound on the evolution speed that 
works for any Hamiltonian, be it Hermitian or non-Hcrmitian. The bounds derived 
here may not be tight for some Hamiltonians and/or for some initial conditions, but 
this statement is equally true for the Hermitian Fleming bound. 

This leads directly to the second goal of this article: to show how to construct 
Hamiltonians for which the evolution speed of the state of interest reaches the 
upper bound for any instant of the evolution. We call these Hamiltonians "maximal 
efficiency" Hamiltonians (or maximally efficient Hamiltonians). For every state 
evolution there exists a family of Hamiltonians that are maximally efficient. This 
family contains both Hermitian and NH Hamiltonians. 

Our third goal is to explore a very special Hamiltonian in this family which is of 
rank-one, non-diagonalizable and similar to a Jordan block with zero-eigenvalue. This 
Hamiltonian corresponds to a NH degeneracy called exceptional point (EP) which 
has only one (geometric) eigenvector [TJ [2TJ [M] • We show that any state evolution 
can be generated solely by such NH degeneracies yielding an EP-driven evolution. 
This special evolution minimizes the Hilbert-Schmidt norm Y_\ ■ of the matrix 

Hamiltonian %. 

We note that the second goal strongly differs conceptually from the so called 
quantum brachistochrone problem. The quantum brachistochrone problem consists 
in finding the (time independent) Hamiltonian which evolves some predefined initial 
state into some predefined final state in a minimal time. This problem was the subject 
of intensive studies during the last few years for Hermitian systems [35} [3BJ [37] as well 
as for NH ones E01 ESI [Ml HOI IIH IH1 H31 HU • As shown in [35JIM1S5] for Hermitian 
systems, the corresponding minimal-passage-time trajectories correspond to geodesies 
in projective Hilbert space. For the evolution problems we are investigating here, the 
trajectories in projective Hilbert space are predefined and not necessarily geodesic. 
Instead, we are searching for Hamiltonians capable of producing evolution processes 
which exactly follow these predefined trajectories with minimal resources. That is, 
we look for efficient evolution and not for a fast evolution. In fact, our optimization 
problem is closer in spirit to the reverse engineering approach used to quicken adiabatic 
evolution H7J HE] ■ Yet there are two main differences: The first difference is that 
in our case we seek only Hamiltonians which yield maximal efficiency. The second 
difference is that we take as input only the evolution of a single state (the state of 
interest), while in [151 ITF1 B5] the number of states needed to be specified is equal to 
the Hilbert space dimension (number of levels in the system) . 

The article is organized as follows: Section [2] contains some basic facts on the 
evolution speed in projective Hilbert space. In section 02 the Hilbert-Schmidt norm 
and the spectral norm of a Hamiltonian are introduced as upper bounds on the 
evolution speed. In section |4j the concept of speed efficiency is introduced, and 
for the predefined evolution of a given state a family of Hamiltonians is constructed 
which ensure a speed efficiency of 100%. The generic properties of maximally efficient 
evolutions are explored. Section [5] is devoted to the special case of a maximally 
efficient evolution which is driven by a Hamiltonian at an EP (a NH degeneracy). In 
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Appendix 1, for completeness we briefly discuss the relation between Bloch sphere 
and projective Hilbert space. In Appendix 2, the norm speed bounds are illustrated 
for a Hamiltonian that describes two optical systems recently studied. In Appendix 
3, the norm speed bound is applied to the matrix Hamiltonian of a VT— symmetric 
quantum brachistochrone. 



2. Preliminaries - the evolution speed in projective Hilbert space P(io) 

Let E S) = C N be a solution of the time-dependent Schrodinger equation (TDSE): 
idt |V) = H(t) |V) , (1) 

where H(t) ^ H*(t) G t£ NxN is the matrix of the corresponding time-dependent NH 
Hamiltonian. Defining the bra-vector (ip\ in a standard wajQas (V| = the adjoint 
TDSE has the form 

-id t {4 ) \ = {^\H\t). (2) 

Our main interest is to study the rate at which states evolve into different states. 
Phase evolution is irrelevant for this purpose. It makes sense, then, to study the 
evolution of states in a space where the phase is eliminated. The so called projective 
Hilbert space (PHS) is exactly suited for this purpose. The well known Bloch sphere 
for two- level systems is closely related to PHS (see Appendix 1), but strictly speaking 
it is not a PHS. For the reader not familiar with PHS we provide a simplified and very 
limited presentation of the basic ideas needed to understand the present work. For a 
more complete and rigorous treatment see, e.g., |49| . 

The angle between two complex vectors |Vi) , IV2) S t C N can be obtained from 
the standard inner product of the two vectors: 

cos e= - «^>l (3) 

The angle acts as a measure of distance between two states: = means the 
two states are identical up to a complex factor and = 7r/2 indicates the states are 
mutually orthogonaQ. 

Now imagine that a state is infinitesimally changed from to |V) + MV) where 
(VI V) 3* (dtp\dijj). The angle between the original state and the modified state, d®, 
can be obtained from Keeping leading orders in \dip) and dO we get: 



d e> = ^tm _ im <y > = dsh . (4) 

(VIV) (VIV) 



ds FS is known as the Fubini-Study metric [33] which describes the length of an 
infinitesimal arc traced on a unit hypersphere by changing a state by \dip). To quantify 
the rate at which states change we will look at the evolution speed defined by: | ^ | 
(or equivalently | ds £ s | ) , which can be interpreted as angular speed/frequency. This 
hypersphere is related to the projective Hilbert space associated with the Fubini- 
Study metric. All states which differ by a complex number are mapped to the same 
point on the hypersphere (hence phase is immaterial in this space). The details of 



1 ) Unlike other choices often made for NH Hamiltonians in order to exploit the bi-orthogonal 
relations of the eigenstates [T] . 

2 ) This is different from the Bloch sphere construction discussed in Appendix 1, where orthogonal 
states correspond to the angle of it between antipodal points on the sphere. 
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this mapping are not important for the present article. What is important is that 
the distance between two states on the hypersphere which differ by \dip) is given 
by ([!]). Formally, the projective Hilbert space of an N-level system is denoted by 
¥(Sj) = CP^- 1 = Cf/C where Cf = C w - {(0,0, . . . ,0)}, C* := C - {0}. As the 
state \ip) evolves in time it traces a certain trajectory in F($j) (i.e. on the hypersphere 
associated with the PHS). We denote the trajectory induced by \ip) by ir{\ip)) G P (■£))■ 
Notice that for any complex function of time, c(t) (c(i) ^ 0): 

7r(|V)) = 7r(c(t) (5) 
Next we wish to establish a relation between the evolution speed | ds Jj. s | and the 
Hamiltonian. Making use of the TDSE (JlJ and its adjoint ©, and introducing the 
normalized state vectors |\&) := \tp) / yj (V'IV')i we find the squared evolution speed in 
¥($)) is given by: 

ds 



dt 



) = (*|Ht( t) ^ (i) |*) _ (*\Hi{t)\*)(*\H(W) =: K{t). (6) 

Henceforth, we refer to K(t) as "kinetic scalar" because it plays a structurally similar 
role like the kinetic energy in classical mechanical systems. The expression © is a 
straight-forward generalization for non-Hermitian Hamiltonians of the corresponding 
evolution speed discussed in |45[ |4U] for Hermitian systems. We note that for those 
systems K (t) just reduces to the instantaneous energy variance K(t) = (\f , |ii" 2 (£)| 1 I r ) — 
(<I'|i7(i)| , I') 2 = AE 2 (t). Further insight into this expression can be obtained by 
introducing an instantaneous orthonormal basis set {\k)}^ =1 G S) = C N , with |\E') 
identified with one of its elements |^) = \j), the kinetic scalar ([6|) in this basis set 
takes the form: 

N 

K = "£{j\Hi\k)(k\H\j) - (j\rf\j){j\H\j) = ^ \(k\H\j)\ 2 > 0. (7) 

k=l k=£j 

Obviously, K is characterizing the total rate for transitions from the given state 
W = \j) to other states of the system. Splitting off the trace of the Hamiltonian 

H = U + fiI, M := Tr(H)/N (8) 

one immediately sees that K is invariant with regard to trace shiftiO, K[H] = K[H}. 
Hence, we can subsequently restrict our attention to traceless Hamiltonians H. In 
Appendix 1, we discuss the dynamics on the PHS-related Bloch sphere and obtain its 
relation to the kinetic scalar. 



3. Upper norm bounds on the evolution speed 

The absolute values of the matrix elements of Hermitian or non-Hermitian 
Hamiltonians are directly defined by the intensities of interactions and the strength 
of the corresponding fields. Naturally, a model remains valid only within the region 
of applicability of the corresponding underlying theory and/or the applicability of the 
approximations made in the derivation of the model. Hence, it is natural to ask for the 
maximal evolution speed achievable by a given quantum system when the resources 
are limited. 

3 ) A time-dependent trace tr(H) = Nfi(t) £ C can be removed from the Hamiltonian H by the 
transformation \ip) — > e ! Jo \tp). Since this transformation involves only a multiplication by 

a complex function the motion in the PHS is not affected by this transformation l(5j- Setups with 
non- vanishing complex traces have been considered, e.g., in |38| . 
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3.1. The Hilbert- Schmidt-norm upper bound on the evolution 

To obtain the first simple upper bound on the evolution speed we notice that: 

K < (^I^KI*) < tr(n f H) = \\U\\% S . (9) 

where ||%||#,s is the Hilbert-Schmidt norm of the Hamiltonian [50]. It is also known 
as Euclidean norm, l^— norm, Schatten 2-norm, Frobenius norm and Schur norm [50J. 
It was possible to use the trace in the last inequality since T-Ol-L is a positive operator. 
Equations ([6]) and © set a bound on the evolution speed: 

^=^<II«|| M . (io) 
Upon writing the HS norm explicitly in terms of matrix elements 

ii^iiL = Ei^i 2 ' ( n ) 

it becomes clear that the evolution speed is limited by the size of the Hamiltonian 
elements and not just by the eigenvalue difference. This becomes very important in 
the vicinity of NH degeneracies as shown in Appendix 3. 



3.2. The spectral norm upper bound — a tighter bound on the evolution 
To get a tighter bound on the evolution speed we use the following inequality: 

N 

K < {m^) = E A *K^}| 2 < max(A fc ) ee \\H\\ 2 SP , (12) 

k=l 

where A& > and \k) are the eigenvalues and eigenstates of the matrix V.^'H, 
T-0l-L\k) = Xk\k). \\H\\sp = \Anax(AA;) is known as the spectral norm of T-L (also 
known as Ky Fan 1-norm [SD]). To understand the second inequality in (fT2")l , notice 
that the states {\k)} constitute a complete orthonormal basis set and ( V I / | V I / ) = 1, so 
that the projection sum satisfies J2k 1(^)1 = !■ Thus : J2k=i x k\(^\k}\ 2 is just a 
weighted average of positive numbers Afc. Such a weighted average is always smaller 
or equal to the largest element. 

Obviously, the Hilbert-Schmidt norm ||%||#s and the spectral norm ||H|||p can 
be represented in terms of the eigenvalues Xk ■ This implies the following useful 
relation: 

||W|||p<||W||^<rank(H)||W|||p. (13) 
Therefore, the spectral norm bound on the evolution speed 

d -^-\ = VK<\\n SP (H) 

is always tighter than the Hilbert-Schmidt norm bound. Simple and useful lower and 
upper bounds on ||%|| S p for H(t) £ C NxN (but not on the evolution speed!) are given 
by: 

max(\n id \) < \\H\\ SP < N maxd^jl). (15) 

Finally, we briefly comment on two extremal cases of traceless 2x2 matrix 
Hamiltonians. 
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• For a NH two-level Hamiltonian % which is similar to a Jordan block with zero- 
eigenvalue % ~ ^2(0) = ( g q J it holds rank(H) = 1 so that 



H ~ J 2 (0) 



SP 



\n\ 



HS ■ 



(16) 



This fact will be important later on in sec. [5l 
• For a Hermitian traceless two-level Hamiltonian with energy separation AE, the 



spectral norm is = \AE\ /2 and we obtain: 

ds FS 



dt 



< \ AE\/2. 



(17) 



As discussed in Appendix 1 the Bloch unit vector, n, is related to the Kinetic 
scalar via: |^| = 2\J~K. Therefore the corresponding upper bound on the 
evolution speed over the Bloch sphere in the Hermitian case is: 



dn 


= 2 


ds FS 


~dt 




dt 



< \AE\ 



(18) 



which is known as Fleming bound [33 . That is, for a two- level Hermitian operator 
the spectral bound coincides with the known Hermitian bound. 

For explicit examples of the speed bound we refer the reader to Appendices 2 and 
3. In Appendix 2, a Hamiltonian that describes certain optical systems is analyzed. 
Appendix 3 studies a PT- symmetric Hamiltonian that was introduced in [53], in 
the context of the 'PT-symmetric brachistochrone problem. In the next section we 
introduce the notion of speed efficiency which quantifies how close the actual motion in 
¥(Sj) is to the speed bound just derived. Later we show how to construct a Hamiltonian 
which reaches the spectral bound for a given motion in ¥(Sj) at all times. 



4. Speed efficiency of quantum evolution 

In this section we introduce the notion of a maximally efficient evolution. We 
wish to compare the actual speed of motion in the projective Hilbert space P(io) 
to the speed bound given by the spectral norm characterizing the available 

resources of the system. We use ||K|| S p since it is tighter than the Hilbert-Schmidt 
norm Let \if>) £ C N be a time-dependent state in an TV-level system that 

induces some predefined evolution tt(| V 7 )) m the corresponding projective Hilbert space 
V(Sj) = CP*- 1 3 n(\ip)). We define the efficiency to be: 

#,W)^<1. (19) 
\\ri\\sP 

It is important to realize that this efficiency is an instantaneous (or local) property of 
H and its solution \ip). The shape of the curve in P(fj) alone has nothing to do with 
efficiency. For example, a geodesic in P(f)) can have efficiency smaller than one, and 
on the other hand, non-geodesic curves can have 100% efficiency. 

Loosely speaking, the value of r/ quantifies to what extent the Hamiltonian 
really uses all its resources to generate motion in P(fj). That is why we call an 
(?7 = 1)— evolution, a "maximally efficient evolution". As an example of inefficient 
evolution, consider a spin in a magnetic field which is not exactly perpendicular to the 
spin direction. The part of the magnetic field which is parallel to the spin is wasted 
as it does not contribute to the precession motion. As we shall demonstrate in the 
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next section, this inefficiency can be fixed by making the Hamiltonian time-dependent 
(rotating the magnetic field in time). In the NH case, reaching 100% efficiency becomes 
even more difficult. As explained at the end of section^ for NH systems the condition 
■j^H = does not guarantee a constant evolution speed, i.e. | L 3 | ^ const. On the 
other hand, the spectral norm is fixed if jj^H = 0. Equation (fH?|) , then, implies that 
rj varies with time and, therefore, the evolution cannot be maximally efficient at all 
times. 

In the next subsection we show how to construct Hamiltonians that are designed 
to generate maximally efficient evolution for a given predefined motion in projective 
Hilbert space at all times. We will demonstrate that such an (77 = 1)— evolution can 
be either Hermitian or non-Hermitian. 



4--1- Maximally efficient evolution 

Our goal in this section is to find a Hamiltonian Hq that generates the same motion 
Tr(\ip)) in P(£j) as but with 100% efficiency. The solution \m) corresponding to 
the maximally efficient Hamiltonian Ho may differ from \ip) only by a time dependent 
complex factor. In short, we look for \m) and Ho that satisfy: 

7r(|m)) = nQii)) (20) 

id t \m) = Ho \m) (21) 

V(t) = 1. (22) 

The first requirement is that the states |m) and \ip) have the same motion in P(Sj). 
The second requirement states that {Ho, \ m )} satisfy the TDSE, whereas the third 
requirement simply means that we are searching for maximal efficiency. To satisfy the 
first requirement we set: 

\m) = c(t) |V>> (23) 

where c(i) is a complex differentiable function of time (see eq. ([5])). For reasons that 
will become clear shortly, we fix c(t) by choosing |m) to be normalized to unity and 
to be parallel transported.]: 

(m\m) = 1, (24) 

(m\d t \m) = {d t (m\)\m) = 0. (25) 

The first condition determines |c(i)|, and the second one yields the phase of c(i) (up to 
a time-independent constant). To find the Hamiltonian that drives |m) with maximal 
efficiency we choose the following ansatz: 

Ho = i \dtm) (m\ — ig \m) (dtm\ , (26) 

where, in general, g can be time-dependent. For g = 1 the Hamiltonian is Hermitian. 
Notice that \dtm) = dt \m) is not normalized and not parallel transported. In 
contrast to |m), \dtm) is also not a solution of the TDSE (with Ho as Hamiltonian). 
However, |m) and \dtm) are mutually orthogonal by virtue of the parallel transport 
we imposed. By applying ([26)) to the state |m) we see that the requirement (|21[) is 
immediately satisfied. To fulfill the remaining third requirement we note that in the 
basis {\m) , \dtm)} the Hamiltonian Ho has only off-diagonal elements so that Ha is 



4 ) If |x) is a normalized state, (xlx) = b then its parallel transported form is given by|x) = 
e — Jo (^l a t' x)dt' ^ |^ satisfies (x 1 9fx) = (dt'X\xJ = 0- This fixes the phase of the state up to a 
constant determined by the choice of the lower limit of the time integral. 
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traceless by construction. To calculate the efficiency we first calculate the spectral 
bound and then the kinetic scalar. Evaluating HqHo we g e t 

h\Ho = (d t m\dt to) |to) (to| + \g\ \d t m) (d t m\ . (27) 

The eigenstates of H^Ho are \ m ) an d \9tTn) and the corresponding eigenvalues are 
{d t m\d t m) and \g\ (d t m\d t m), respectively. The spectral norm is given by: 



\\Ho\\sp - V(dtm\d t m) max(l, | 5 |). (28) 
Using the fact that (n \Ho\ n) — and Eq. © and (f2"7| we get that: 

K (H Q , |m» = {d t m\d t m) . (29) 
The efficiency, then, is given by: 



= y/K{U ,\m)) = V(d t m\d t m) 

WHoWsp \J (d t m\d t m) max(f, \g\) 

Clearly, maximal efficiency r/ = 1 (third requirement (|22D ) is achieved provided that: 

lfll<l- (31) 
In summary, given any arbitrary state equations (|2"3"|) - (f2"6|) together with (f3"Tjl 

show how to construct maximal efficiency Hamiltonians that generate the same motion 

in projective Hilbert space P(Jj) as 

It is instructive to look on the instantaneous eigenvalues of Hq . From 

Ul = g \d t m) (d t m\ + g (d t m\d t m) \m) (m\ (32) 

and tr('Ho) = it follows that the instantaneous eigenvalues E±(t) of Hq are: 



E ± {t) = ±y/gy/(dtm\dtm). (33) 

In case of g £ R, these eigenvalues are real for g > and purely imaginary for g < 0, 
i.e. E±(t) £ RU iM.. This indicates a hidden instantaneous pseudo-Hermiticity of Ho 
— in a certain analogy to the consideration^ in Appendix 3. One of the key points 
of this work is that the evolution can be maximally efficient regardless of whether the 
Hamiltonian is Hermitian or not. 

Finally, we note that the Hamiltonian Ho in shows some structural analogy 
to the brachistochrone Hamiltonians for Hermitian systems (constructed in [36 ). 
In fact, Ho extends the geodesic-trajectory paradigm of |35[ [36l [45] to maximally 
efficient evolution regimes over arbitrarily predefined time-dependent trajectories in 
P(fj). Moreover, the constraint AE — const is replaced by the constraint i] = 1. 



4-2. Inherent properties of the maximally efficient evolution 

Here we wish to highlight three points which are generic for maximally efficient 
evolutions. The first point concerns the fact that \m) is normalized and parallel 
transported. In the construction of Ho we demanded (m(t)\m(t)) — 1 and (TO|(9 f m) 
(). Now we wish to show that if these constraints are relaxed the spectral norm will 
increase even though the motion in P(fj) remains unaltered. According to (fT!?|) the 
efficiency will drop below 100% by this modification. Assume we wish to change the 
amplitude and phase of (to) by some complex factor e _4V '*' where (p(t) is some complex 
number. This is accomplished by adding Ho a diagonal term so that: 

H new =H + \m) (m\ d t ip(t). (34) 



5 ) We leave a corresponding detailed investigation to future research. 
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To keep the trace zero, another diagonal term must be added as well, in principle, but 
it is of no importance to the present discussion. This transformation does not change 
the value of K . The spectral norm squared is the largest expectation of r hi\ lew l-i new , 
so: 

> ( m \^-ne W ^-new\ m) = \dtf(t)\ 2 + (d t m\d t m) > (d t m\d t m) (35) 

Since K remained the same and spectral norm increased, we see that the efficiency is 
now: 

y/(d t m\d t m) y/(d t m\d t m) 

Vnew = -rr^ |j < , < 1, (36) 

ll«neu,|| SP ^\d t( p(t)\ 2 + (d t m\d t m) 

This decrease in efficiency with respect to Ho expresses the simple fact that changes 
in phase and/or amplitude also require spectral norm resources from the Hamiltonian. 
In order to direct all resources to motion in P(io), any phase and amplitude changes 
should be avoided. 

The second point concerns the role of g. While the Ho found earlier conserves 
the norm of |m), it doesn't do so for other initial states (with the exception of the 
Hermitian case g = 1). Moreover, while \m) evolves exactly in the same way for all 
values of g, the evolution of other states strongly depends on the value of g. This is 
demonstrated in the example shown in Figure 1. The state of interest in this example 
was chosen to be: |m) = cos(t) |f) +sin(t) We considered two different maximally 
efficient Hamiltonians. The first one is Hermitian (g=l) and the other is not (g=- 
0.8). The Hamiltonian is constructed using the recipe in Sec. 14. ft If the initial 
state is \ip(t = 0)) = \m(t = 0) = |t) we observe that, as expected, both Hamiltonians 
generate the same evolution. Yet, if \ip(t = 0)) = \\), the different Hamiltonians 
generate different evolutions and the effect of g becomes apparent. 

The third point that we note is that for periodic motion in ¥(£)) the state 
| to) accumulates only an Anandan-Aharonov phase |5Ij . since the dynamical phase 
(to \H\ to) is zero for maximally efficient evolution. Once again, this is the result of 
wasting no resources on phase accumulation. 



5. EP driven evolution 

A special case of great interest is g = 0. Equation ([32]) shows that in this case H\ = 0. 
This can only happen if Ho is similar to a rank-one Jordan block with zeros on the 
diagonal, i.e. when Ho ~ ^(O). That is, Ho(g — 0) describes an EP operator — a NH 
degeneracy. The dynamics can be fast even though at each instant the instantaneous 
eigenvalue difference is zero. This degeneracy is time- dependent. Moreover, the 
orientation of the single geometric eigenvector of a Hamiltonian Hq ~ ^2(0), associates 
a preferred directionality to this degeneracy. The directionality of the EP at each 
instant of time is chosen such that it induces the desired dynamics. Thus it appears 
natural to name this evolution an EP Driven Evolution (EP-DE). Another unique 
feature of this evolution can be seen by evaluating the Hilbert-Schmidt norm. For a 
general value of g the HS norm is: 

\\Ho\\ 2 HS = (l + \g\ 2 )(d t m\d t m). (37) 

Obviously, the HS norm takes its minimal value for EP-DE (g — 0), i.e. from all 
the possible maximally efficient evolutions the EP-DE has the minimal HS norm. 
At this point the HS norm ||Wo||jjs and the spectral norm ||"Ho||sp coincide, a fact 
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Figure 1. (color online) The value of g has no effect on the evolution when 
applied to the initial state "Ho was designed to propagate efficiently (North pole 
of the Bloch sphere in this example). Yet when applying Wo to a different initial 
state (South pole) the g may completely change the evolution. The large black 
dots mark the starting points of the evolution. See text for details. 



mentioned in (IT51) . Equation (fTTj) shows that the EP-DE provides the minimal value 
of \^-ij\ 2 f° r a gi ven trajectory in P(ij). 

Conclusion 

The concept of speed efficiency was defined using spectral speed bounds derived for NH 
or Hermitian systems. A recipe for the construction of 100% efficiency Hamiltonians 
for any evolution in the projective Hilbert space was given. These Hamiltonians 
contain a free parameter, "<?". 100% efficiency is obtained for \g\ < 1. The Hermitian 
case corresponds to g = 1. We conclude that it is possible to have a Hamiltonian which 
is both 100% efficient and Hermitian. The g — case corresponds to a Hamiltonian 
which is not diagonalizable, i.e. the evolution is driven solely by a time-dependent 
NH degeneracy (exceptional point). This particular evolution minimizes the quantity 
Si j \T~t-ij\ 2 with respect to all other 100% efficiency evolutions considered in this work. 

Appendix 1 - The Bloch Sphere and the Fubini-Study metric 

The dynamics of NH 2x2 matrix systems in f) = C 2 9 \ip) is conveniently analyzed 
as dynamics on the Bloch sphere. The latter is spanned by the unit vectors 

n {t) = M = e S 2 C M 3 . (38) 
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Its close relationship to the projective Hilbert space CP 1 = P(f)) = C^/C* ~ S 2 can 
be seen by the explicit comparison of the Bloch sphere metric with the Fubini-Study 
metric of CP 1 . For a qubit |W) e C 2 parametrized as |*) = (cos(0/2), e^ sin(6/2)) T , 
6 £ [0,7r], 4> E [0,27r] it holds n = (sin(0) cos (</>), sin(0) sm((f), cos(6»)) T and the Bloch 
sphere metric reads 

dh 2 = d8 2 + sin 2 (9)d(f> 2 . (39) 

For the same state |^/) the Fubini-Study metric ((6]) reduces to 

ds FS = \( d ° 2 + sin 2 (9)d<p 2 ) (40) 

and therefore: 




This is closely related to the fact that orthogonal states are antipodal on the 
Bloch sphere having a geodesic distance n, whereas the corresponding Fubini-Study 
distance as discussed in section 2iss_Fs = = 7r/2. 

The main results of this article can be expressed using the Bloch Sphere and the 
NH Bloch Equations (see for example [5]) formalism, but we found that the results 
are more neatly described by the "ket-bra" operator formalism and the Schrodinger 
equation. Moreover, unlike the NH Bloch Equation formalism the "ket-bra" formalism 
is applicable to a multilevel system without any alterations. 



Appendix 2 - Speed bounds in optical systems 

Let us examine the NH evolution in optical systems where the Hamiltonians are 
explicitly known and a two-level description is either a good approximation or even 
exact. Consider the Hamiltonian Hq introduced and studied in |31) in the context of 
"EP cycling" [TH 

0' (42) 

where z, the propagation coordinate, plays the role of time. This Hamiltonian can 
describe different physical systems. In [TB], it was used to describe the evolution of the 
transverse electric field and its spatial derivative in a waveguide. In this system q(z) 
is proportional to the change in the index of refraction with respect to vacuum. In [5] 
H(z) describes the evolution of the two optical polarizations in crystals. q(z) in this 
case is related to the change in the transverse part of the reciprocal dielectric tensor. 
The eigenvalue difference of H(z) is AE = 2y/q(z), and at q — a NH degeneracy 
(an EP) forms which is experimentally well accessible in these systems. 

From the structure of 'Ho(z) in (|42|) . it is obvious that for q = and AE(q = 0) = 
the evolution speed does not vanish for all states. Using (|4Tjl and the kinetic scalar 
definition ((6j), one finds that the non- vanishing angular velocity is |g | = |^| = 2 \ q(z)\ 
for the spin-up state |f) = (1, 0) T , and \6\ = 2 for the spin-down state = (0, 1) T . In 
particular, at the degeneracy q = the spin- up state becomes an eigenstate (6 = 0), 
whereas the spin-down state still has a speed 6 = 2 regardless of AE = 0. For q ^ 
the evolution speed of a general state (not necessarily spin-up or spin-down) can be 
shown to be limited by 

\6\ <2max(l,|g(z)|) <2^jl + \q\\ (43) 
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where the last inequality simply follows from max(|a I , |b|) < \J\a\ 2 + \b\ 2 . From the 

norms of the Hamiltonian (02) , 1 1 "H 1 1 SP = max(|i| , |— iq(z)\) , \\H\\ HS = \Jl + \q\ 2 , we 
see that 

\6\<2\\H\\ SP <2\\H\\ HS . (44) 

Hence, the maximal speed exactly fits within the spectral norm bound of H(z). Other 
optical systems whose evolution speeds near EPs can easily be studied are discussed, 
e.g., in [15] , 



Appendix 3 - Spectral speed bounds for pseudo-Hermitian and 
VT— symmetric two- level Hamiltonians 

We start from a general type NH traceless Hamiltonian, H, written in terms of the 
Pauli matrices a: 

H(t) = [d(t) + 0{t)] ■ a , (45) 
where a, j3 £ R 3 are some time-dependent real vectors. The eigenvalues are: 

E± =±\J(a + i0) 2 = ±\fd 2 - 2 + i2a ■ 0, (46) 

and become real or pairwise complex conjugate for 

a- 0=0 => E± E R U iR. (47) 

Operators and matrices with this specific spectral behavior are known to be symmetric 
under an anti- unitary transformation [52] , to be pseudo-Hermitian [53] and self-adjoint 
in a Pontryagin space |54| (a finite-dimension type version of a Krein space [55 ). 

For Hamiltonians % with a ■ j3 = the SP norm and the eigenvalue difference 
reduce to 

\\ j H\\sp = \*\ + \0\, AE = 2^a 2 -0 2 . (48) 
Obviously it is possible to have an arbitrary large ||H|| S p and a vanishing energy 
difference by choosing \d\ — > \/3\. This choice corresponds to an EP limit for which 
\AE\ I ||H|| S p — > 0, since the eigenvalues are very small near the degeneracy while 
sp remains roughly constant and finite. 

A simple example for a NH Hamiltonian with a similar type of behavior is 
the Hamiltonian % used in studies of the VT— symmetric quantum brachistochrone 
problem [29] 

. . / zrsinv s \ , . . 

n=[ A . . = sa x + ir sin X <r z , r,s,xeR. 49 

This complex symmetric Hamiltonian is VT— symmetric, [VT, Ti\ = and V— pseudo- 
Hermitian, VT~L = T-OV, with the parity operation given as V = <j x and the time 
reversal, T, as complex conjugation. In |29| . it was shown that for certain parameter 
combinations r, s, x such Hamiltonians with fixed and purely real eigenvalues difference 
AE e R can evolve a given initial state 1^/) 6 C 2 into an orthogonal final state 
\^f) € C 2 , (^iI^f) = 0, in an arbitrarily short time interval. Due to the finite 
geodesic distance ir between these antipodal states l^i) and \ ^f) on the Bloch sphere, 
the corresponding evolution speed should diverge in this limit. The concrete relations 
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can be easily obtained in terms of the simplifying reparametrization r sin \ 
which yields 



H = s 



1 



1 

-i sin a 



AE = 2s cos a, 



\n\ 



SP 



= \s\(l 



sin a 



s sin a 



(50) 



As demonstrated in |39) . the ultra- fast evolution regime predicted in |29) corresponds 
to an EP-limit a — > ±ir/2 so that for fixed AE = const it holds s = AE/ (2 cos a) — > oo 
and 



H 



i 1 
1 -i 



\\n\ 



SP 



2\s\ 



\AE\/\\U\ 



SP 



0. (51) 



According to (|4"Tj) , this would indeed allow for diverging evolution speeds on the Bloch 
sphere 

dn 



dt 



2Vk<2\\H\ 



sp 



oo. 



(52) 



From this diverging spectral norm one might be led to the conclusion that actually 
such ultra-high evolution speeds and corresponding ultra-short evolution times might 
be forbidden by the limited resources of the system and the validity region of the 
model used. Both would set some natural upper bounds (ultra-violet cut-offs) on 
the evolution speed. This would be true if one were keeping within the present NH 
setups. Nevertheless, the same ultra-high-speed evolution regimes can be induced 
in subsystems of entangled Hermitian systems in larger Hilbert spaces [43]. Due 
to geometric contraction effects the corresponding evolution speed of the associated 
(Naimark-dilated) Hermitian system in the larger Hilbert space will remain finite, 
well-behaved and much below any ultra-violet cutoffs. 

For completeness we note that the present evolution speed considerations are 
closely related to questions for possible lower bounds on evolution times (quantum 
brachistochrone problems) and possible violations of such bounds. Corresponding 
intensive theoretical studies for Hermitian setups (35] [36l [37] in the early 2000 's 
have been followed by investigations of various aspects of non-Hermitian systems 
(VT- symmetric [IH [501 131 SEE S3 SI S3] and quasi-Hermitian ones [H], as well 
as other of more general non-Hermitian types |38|h 

VT— symmetric setups (56] [57] [58] have been experimentally studied via special 
arrangements of gain-loss components (active VT— symmetry) and components of dif- 
ferent loss (passive VT— symmetry) in optical waveguide systems [FJHI], microwave 
billiards [TF], electronic LRC-circuits [Ti?l |2"U] and in mechanical systems of coupled 
pendulums [22]. 
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